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1. Introduction and summary of known results 

In this work we consider Supersymmetric Yang-Mills Quantum Mechan- 
ics (SYMQM) [1, 2], which in D = 10 dimensions and for SU(N -► oo) 
gauge group according to the remarkable hypothesis [3, 4, 5] is equivalent to 
M— theory of DO branes. Despite the fact that the three loop calculations 
[6, 5] question the exact equivalence, the SYMQM models possess a lot of 
fascinating properties required in M— theory, i.e. continuous spectrum of 
scattering states with the threshold states describing a supergravition [7]. 

The SYMQM emerges from the dimensional reduction of supersymmet- 
ric Yang-Mills field theory in the D — 1 dimensional space to the effective 
Quantum Mechanics of zero momentum modes in a single point [2, 8]. In 
a gluinoless sector the energy spectrum of these system is identical to the 
spectrum of 0— volume glueballs [9, 10, 11, 12]. Although the resulting mod- 
els are much simpler than the original field theories, they are rather complex 
with non-trivial solutions. Researches on the systems with D < 10 for var- 
ious gauge groups provide the global understanding of whole family. Thus, 
the SYMQM provides a simple laboratory to study many properties of su- 
persymmetric systems [2, 8]. For D = 2 and SU(2) gauge group the system 
is exactly solvable. Going further to D = 4 dimensions the model becomes 
non-trivial. It possesses both localized and non-localized states. 

The simplest method of solving the energy eigenproblem for the D = 4 
model with SU(2) group was introduced in Refs. [12, 13]. According to 
this method acting with bosonic and fermionic creation operators on the 
empty state one can construct a basis of Fock space. Next, in order to find 
the spectrum and the eigenfunctions of the Hamiltonian one can solve the 
eigenproblem of the Hamiltonian matrix calculated in the constructed ba- 
sis. Since we have an infinite number of the states in the Fock space in 
order to perform numerical calculation we have to cut some basis vectors 
off. Therefore, in the above method one considers only Fock states with 
a number of bosonic quanta tib smaller or equal to cut-off B. At the end 
making analysis of the spectrum as a function of cut-off B we are able to find 
the energy for the model without a cut-off. It turns out that due to com- 
mutativity properties between the Hamiltonian and the operator of number 
of fermionic quanta, hf, the system separates into sectors with determined 
np = 0, 1, . . . , 6. Thus, we can solve the eigenproblem in smaller Fock ba- 
sis with fixed np. The above method allowed to calculate the spectrum in 
all fermionic sectors with an arbitrary angular momentum j for the cut-off 
B = 8 and the Witten index for the model [12, 13]. 

It turns out that the Hamiltonian also commutes with an operator of 
the angular momentum j. This makes possible to construct the basis with 
fixed both np and the angular momentum, j. It was done by van Baal in 
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Refs. [14, 15] for n F = 0,2,4,6 and j = but with the cut-off B = 39. 
Calculations for other sectors are more complicated within this approach. 

The aim of this work is to analyse the spectrum of the D = 4 model with 
SU(2) gauge group more precisely by enlarging cut-off yet further. Using 
the van Baal's method we concentrate on the sector with F = 2 and j = 
where the supersymmetric vacuum state appears to be [15]. In this sector the 
eigenvalues of the Hamiltonian form discrete as well as continuous spectrum 
[16, 17]. For high cut-offs the wave-functions of discrete states converge with 
B and therefore we are able to build these eigenstates and describe properties 
of their wave- functions. Moreover, we can test the scaling properties [18] 
appearing for the continuous spectrum. 

In the beginning of the next section we introduce a notation following 
Ref. [15] and describe the SYMQM model and show how one can construct 
Fock space in the F = 2 an d j = sector and calculate the Hamiltonian 
matrix. We also present the construction of the eigenstates of the Hamilto- 
nian. Next we show the results which include the energy spectrum of the 
model for cut-off B < 60. Moreover, we construct and describe properties 
of the eigenfunctions for the discrete and continuous spectrum. At the end 
we confirm the scaling relation for the continuous spectrum which was de- 
rived in Ref. [18] for the cut free-particle and only conjectured for interaction 
systems. 

2. Hamiltonian eigenproblem 

2.1. Definitions 

Following Refs. [9, 15] the supersymmetric Hamiltonian of Yang-Mills 
quantum mechanics with the SU(2) gauge group has a form 

H6 aft = ±{Q a ,Ql}, (2.1) 

where Q and are the generators of SUSY in the coordinate representation 
of the field theory defined as 

Qa = (-^ - iB?) Q a = + > (2-2) 

with Pauli matrices o~i = r J and Weyl spinors Aa- Here, a colour magnetic 
field of SU{2) is defined as 



(2.3) 
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where V? are bosonic variables with colour indices c = 1,2,3 and spatial 

indices j = 1,2,3. The Weyl spinors Xa with j3 = 1,2 spinor indices and 
a = 1,2,3 colour indices satisfy the following anti-commutation relations 

{\ aa ,\ b P} = a% a 5 ab , {X aa ,X b/3 } = 0, {A a ",A^} = 0, (2.4) 

where a is the unit matrix. 
The SUSY generators satisfy 

{Q a , Qa} = 2(a ) aA n - 2{o i ) a6 y l a g a , (2.5) 

where 

is the generator of the SU{2) gauge transformation while the Hamiltonian 
density [8] is defined by 

W = -I , T ,L n + IzfB? - ige abc \ a a j \ b V^ . (2.7) 
2 dVf-dVf- 2 1 1 3 y ' 

One may rescale bosonic variables V® introducing new variables cf: 

Thus, performing the approximation of constant fields with g(L — ► 0) — > 
[9] with g = g(L) L 3 we obtain a Hamiltonian H independent of the coupling 
constant g: 

' d 3 xH = g 2/3 (L)H/L, (2.9) 



where 

H = H B + H F , (2.10) 

with the bosonic part 
and the fermionic one 

H F = -ie abd \ a a l X b cf. (2.12) 
Here, the colour magnetic field is rescaled to 

Bf = --£ijk£abdc b jcf ■ (2-13) 
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2.2. A straightforward approach 

In this model there are 6 fermionic degrees of freedom, Af . An operator 
which gives the number of fermions, np, commutes with the Hamiltonian 
(2.10). Thus, our system splits into 7 sectors enumerated by np values. 
Since the sectors are related by a particle-hole symmetry 

rip <-> 6 — rip , (2-14) 

we have four independent sectors described by np = 0, 1, 2, 3. 

In order to solve the eigenequation of the Hamiltonian (2.10) we con- 
struct an infinite basis of the Fock space by acting with cf and on the 
empty state: 

l»>= E ^...c-At +1 ...AfjO>, (2.15) 

contractions 
{ai ,a r } 

where df and A£ are given by appropriate linear combinations of the creation 
and annihilation bosonic and fermionic operators, respectively, while the 
sum goes over gauge invariant linear combinations of cf and A£, [12] 1 . Next, 
we act with (2.10) on the basis states (2.15) calculating matrix elements of 
the Hamiltonian. Finally, we diagonalize the Hamiltonian matrix finding 
eigenvalues and eigenstates of (2.10). 

Contrary to fermions we have an infinite number of bosons, tlb = 0, . . . , oo. 
Moreover, the operator which describes a number of bosons tib does not 
commute with the Hamiltonian (2.10). Therefore, our Hamiltonian matrix 
is infinite. To simplify the problem we cut the basis (2.15) by considering 
only the states with hb < B, where B is defined as a cut-off. Next, analysing 
the limit B — > oo we recover the spectrum of the full, infinite Hamiltonian 
matrix. 

Eigenstates of the model are linear combinations of basis states (2.15). 
It turns out that for given hb we have a large number of states. Even if we 
divide the problem into the separate fermionic sectors, the separate bases 
have large numbers of states. One can see it in Table 1 from Ref. [13] where 
the authors show the size of the basis containing states generated in specified 
fermionic sectors labelled by np and for given n#. In this Table N s describes 
a number of basis vectors with a given bosonic quantum numbers whereas S 
is a sum of vectors in the basis with maximal quantum number tib- The last 
column shows that the number of bosonic states S# = So + S2 + S4 + Sq 
is equal to the number of fermionic states S^ = Si + S3 + S5 as required 
by supersymmetry. 

1 Here, calculations are performed in Fock-space representation 
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5 
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445 


240 540 
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91 
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600 1140 
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112 


336 672 
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1490 


720 1860 





8 


126 


238 


336 1008 
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2450 


1500 3360 





Jmax 




8 


17/2 




9 


19/2 





Table 1. Size of the basis generated in specified fermionic sectors np from Ref. [13] 



E 




345678 



345678 



2 4 6 8 



2 4 6 8 



riF 



= 



1lF 



= 1 



=2 



riF 



=3 



Figure 1. Dependence of the energy £onaa cut-off B > ns in sectors with various 
tif obtained in Ref. [13] 



In Fig. 1 we present a dependence of the energy on the cut-off. One can 
see that for this large number of bosonic quanta the energy converges to 
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the constant value rather slowly, especially for np = 2,3,4. Thus, in order 
to simplify this problem, we find another operator which commutes with 
the Hamiltonian, i.e. the total angular momentum j. Using the operator 
describing a fermion number, Uf, and the total angular momentum, we can 
divide the spectrum into smaller sectors labelled by np and j. This will be 
discussed below. 

2.3. Analytic change of variables 

In the rest of the paper we consider the case with np = 2 and j = 0. 
This case is most interesting because it contains the vacuum state of the 
system. Moreover, the energy has discrete and continuous spectrum. We 
will use the notation following the work of van Baal [15, 14]. 

The gauge invariance and the vanishing angular momentum j allow us to 
reduce a number of bosonic variables [15]. The reduction can be rewritten 
as a diagonalization of cf : 

3 

<■■; Y, n >j-/i' j "- (2-i6) 

3=1 

where matrices R,T £ 50(3). Thus, from nine bosonic variables c" one 
obtains three invariants of both gauge and rotation groups: 

— oo < Xj < oo , where z = 1,2,3. (2-17) 

Additionally, the second set of invariant variables (r, u, v) is introduced 
where 

r 2 = (c?) 2 = , (2-18) 

3=1 

defines the radius r in the {xj}— space, so that < r < oo, while 

3 

i>j=l 

takes values < u < | and corresponds to the rescaled bosonic potential 

V(x) = x\x\ + x\x\ + x\x\ = ur A . (2.20) 

An example of equipotential surface is shown in Fig. 2. The potential has 
minimum on x,— axes. Due to the rescaling property (2.20) of V(x) all 
equipotential surfaces with different values of energy have the same shape. 
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Figure 2. An example of equipotential surface for V{x) — ur 4 — x\x\ + x\x\ + 
x\x\ = 0.1 



It turns out that for np = 2 the fermionic part of the Hamiltonian breaks 
the symmetry 

X1X2X3 —X1X2X3 . (2-21) 

Thus, in order to describe the system completely one has to introduce the 
remaining variable 

3 

v = r' 3 det c = r~ 3 J| Xj , (2.22) 
3=1 

which may take positive as well as negative values —l/y/27 <v< l/y/27. 

One can notice that the transition from (r,u,v) to (xi,X2,x^) is not 
unique. The transformation functions Xj(r, u, v) between coordinates are 24- 
valued functions and we can choose one branch, i.e. imposing |xi| < x% < X3. 
It turns out that for uf = 2 and j = bosonic terms of our calculations can 
be rewritten in (xi,X2,a;3) variables or equivalently in (r,u,v) coordinates. 
After integrating out six angular variables, the integral measure changes into 



d 9 c = ^TT 4 Jd 3 x , 



(2.23) 
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with the Jacobian J defined as 

J = II \ %2 i - x )\ = r V^ 2 (! - 4 «) - v 2 {4: - 18u + 27v 2 ) . (2.24) 

i>j 

Let us consider the basis vectors (2.15). For tif = 2 they contain two 
fermionic creation operators A^. We can form either antisymmetric 

|V) = Vj c l j c = -IzV^SatcXU^f^lO) , (2.25) 

or symmetric 

\S) = S ab J ab = -S ab \%\yP\0) , (2.26) 

combinations, where Vj C and S ab are antisymmetric and symmetric com- 
binations of bosonic variables c" respectively, while is related to Pauli 
matrices by = \rj. The covariance gives us their structure as 

V/ = Mr, u > v Y c V r ~ h ^ u > v)B«/r 2 + h 3 {r, u, v)c b c b k c a k /r 3 , (2.27) 

and 

S ab = fu(r, u, v)5 ab - h 5 (r, u, v)c^c b /r 2 + h 6 (r, u, v)c b c d c d k c b k /r 4 , (2.28) 

where hi(r, u, v) are functions of invariant variables. Finally, our basis vector 
may have both invariant parts, so that 

|*> = |V> + |5>. (2.29) 

The Eqs. (2.27)-(2.29) can be rewritten as 

6 

\^) = J2h li {r,u,v)\e ll (u,v)), (2.30) 

where 

|ei («,«)) = cyrl\, \e 4 (u,v)) = 5 ab J ab , 

\e 2 (u,v)} = B«/r 2 l\, \e 5 (u,v)) = c«c b /r 2 J ab , (2.31) 

\e 3 (u, v)) = c b c b c a k /r 3 l\ , \e 6 (u, v)) = c b c d c d k c b k /r 4 J ab . 

Lowering of indices can be done with e a p = = — IT2 for spinor indices, 
S ab for colour group indices and 5ij for space indices. Now, we are ready to 
reformulate our Hamiltonian (2.10) in terms of (r,u,v). 
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The fermionic part of Hp can be expanded in the space spanned by 
M = |e M (u,v)), i.e. 



(2.32) 



which gives a 6 x 6 matrix in a form 



( 
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— V 
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1-u 
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1 








—v 







-1 








-1 


-1 






2 


Av 


2-4ti 















-2 



























-2 











/ 



(2.33) 



The vectors \e^(u,v)) are not orthogonal and their scalar products 2 can be 
written in a matrix form as 



N^ = (e,\e u ) F 



( 


8 


2Av 


8X 








° \ 


24t; 


8u 


8v 













8X 


8v 


8Y 






















12 


4 


AX 













4 


AX 


AY 


V 











AX 


AY 


AZ ) 



(2.34) 



where 

X = l-2u, Y = l-3u + 3v 2 and Z = 1 + 2u 2 - Au + Av 2 . (2.35) 

Thus, calculating scalar products of the states (2.30) we have to use the N up 
matrix, i.e. 



f d 9 c Y, Kn^k, 

J n, V =l 



(2.36) 



with h v defined in (2.27)- (2.28). Similarly, the matrix elements have a form 

6 

= / d"c y h* H up h! . (2.37) 



<P\H\V) = f<fc £ K 

H,u,p=l 



The kinetic part of the Hamiltonian (2.10) is given by 

i d 2 



2 {del) 2 



1 7-1/ ^ 9 Tf-\ 9 1 

~2 J {x) d^ J{x) d^ = -2 [! 



dr dr r 2 



(2.36 



2 These scalar products, i.e. are performed with integration only over fermionic 

degrees of freedom 
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where the Laplacian on a 9— dimensional sphere reads 

q2 q2 q2 

A(u,v) = 4(3v 2 + u-4u 2 )-— 2+8(1 -3u)v—— + (u-9v 2 ^ 



(du) 2 dudv (dv) 2 

+4(2- llu)— -30u— . (2.39) 
ou Ov 

Similarly to (2.32), the matrix elements of A(u,v) can be calculated as 

^6 6 

u, v)\e^ (u,v))= \eAu,v))H^h^r,u,v), (2.40) 
where H 1 ^ is found in Ref. [15] as 



with 



(2 — 4u)d u — 3vd v d v + 2vd u 3vd v + 6ud u 

A v = 2 [ d v (2 - 8u)d u - 6vd v -6vd u 

—2d u —d v —12ud u — 9vd v 

(2.42) 

2vd v -8v 2 d u 

A 1 S = 2\ (4 - 8u)d u - 6vd v Avd v + 8ud u \ . (2.43) 
-4d u -16ud u -12vd v 



and 



-4 7 \ / 3 1 

A v = 2 | 0-9 , A£ = -9 11 j (2.44) 
-15 / \ -22 



To sum up, we obtain the 6x6 Hamiltonian matrix in a form 

= _I,5A«v-8 |_ r 8 1_ + r -2^u + I^ r 4 u + r ^jf (2.45) 
2 or or 2 



where iJ^f = H 1 ^ /r. 



12 



SUSYWFA PRINTED ON 1ST FEBRUARY 2008 



2.4- The Fock space 

Following Ref. [15] we choose the eigenvectors of the harmonic oscillator 
as the basis vectors (2.30). Substituting (r, u, v) variables these eigenfunc- 
tions separate into a spherical and radial part as follows 3 

6 

\¥ n ^(r,u,v)) =^h ( ? Am \r,u,v)\e ll }, (2.46) 

(U=l 

with 

h^ m \r,u,v) = y^(u,v)4> £ n (r) . (2.47) 
Next, to ortonormalize the basis (2.46) the Gram-Schmidt process 

\¥ n ^(r,u,v)) ° rt ^> rm - \¥ n ' e ' m \r,u,v)) , (2.48) 

is performed. 

Construction of the spherical harmonics 

yVri = ( u ,v\e,m) = (# m) ,# m) , . . . ,3f' m) ) , (2.49) 

is shown in Ref. [15]. They satisfy the eigenequation 

H A yV> m ) = L(2L + 7)y^ m) , (2.50) 

where the angular momentum L = — 3) is a half-integer number and its 
degeneration is described for even 2L by 

m = 0,1,2,... ,L, (2.51) 

and for odd 2L by 

m = i,|...,L. (2.52) 

A first few spherical harmonics are presented in Table 2. As we can see the 
spherical harmonics are six-component vectors of polynomials in u and v. 

The radial part of the Schrodinger equation for the harmonic oscillator 
gives the eigenequation 



V^ r+ i (<+4) f- 3) + V 

2 2 r 2 2 



4>l{r) = &J n {r). (2.53) 



3 In the following sections all sums will be written directly. There is no summation 
over repeating indices assumed a priori 
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i(*-3) 



rn 



y(l,m) _ (3;^'™) ^> m ) 3^' m ^ 



0,0,0,1,0,0) 735/2/87T 2 



1,0,0,0,0,0) V105/16tt 2 




1 



0,1,0,0,0,0) 71155/2/16^2 
0,0,0,-i,l,0)3V77/167r 2 



0,0,0,^,0,0)^15015/16^ 
-jl,0, 1, 0, 0, 0) liy27375/32^ 2 



2 
2 
2 




1 

2 



,-^,0,0,0,0) 39\Z77/32tt 2 



0,0,0 10 



nr ^,l)42V7/86/16vr 2 
0,0,0,-^ + n, 



13' 



i,i)3y6T4972/167r 2 



0,0,0,-%,v, 



195' 15' 



0) 3A/510 51/2/ 167T 2 
0, 0, 0) 39VT7^5/64^ 2 

13d 



44 ' 44 



,0,0,0) 33V221/7/167T 2 



3 
3 
3 
3 




1 
2 
3 



10v 
17 



, u, 0, 0, 0) 51^4389/5/32^ 



~TS> -M + u ' lr> °> °' °) 3971 7765/7 /64tt 2 
°> °> °> MS ~ ff + 1,2 > °> °) 663^209/7/32^ 
nnnA_3«_A + ? , 

U ) U ) U ) 51 1 "7 1 I - 



17 



, ^)5iy27T777/327r 2 



0, 0, 0, §§, -*$,v) 969^231/10/32^ 
0,0,0,-^ +uv, 



15v 



10 u I ,,2 
969 19 ^ 17 ' 19 

2w 2 2w 



6_u 38u ' 
13' 65 
2t> 2 



39 v / 53295/27 327r 2 



44 
399 



93m 
133 



+ 



§, -3I3, 0, 0, 0) 969V429/14/327T 2 

^ + it, 0, 0, 0) 19A/51051/2 /647T 2 



7 ' 7 ' 



Table 2. Orthonormal spherical harmonics for L < 4 from Ref. [15] 



The solution to Eq. (2.53) has a form 



t (r) = V2n! 



e- r2 / 2 /- 3 4 + ^(r 2 ) 

r(n + £ + §) 



(2.54) 



where L^(x) are Laguerre polynomials, while quantum numbers 

n = 0,1, 2,..., 00, (2.55) 

and 

£ = 2L + 3 = 3,4,5, ...,00, (2.56) 
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enumerate different radial solutions. The oscillator energy defined by (2.53) 
is a function of the above quantum numbers: 

E e n =e + 2n+l (2.57) 

Since the number of degrees of freedom in the system is nine a number of 
bosonic quanta for a given basis vector (2.48) reads 

n B = E e n -^ = £ + 2n-3. (2.58) 

The above formula will be used to relate the quantum numbers i and n to 
the cut-off B > riB- 

Let us consider the integrals over (r, u, v) which appear in the scalar 
products of the wavefunctions, (2.36) and (2.37). Due to the simplicity of 
(r,u,v) variables and the polynomial form of the eigenstates y n to (2.50), 
one can perform the integrals of the spherical harmonics over u and v: 

X id = [ d 9 c u { v 2j , (2.59) 



Jr=l 

by making use of the recurrence relation 

= 4i(l + i + 4j)X,_ lj + 12i(i - l)Xj-2j+i + 2j(2j - IJXj+ij-i 
hJ (M + 6j)(Ai + 6j + 7) 

(2.60) 

where ^0,0 = 327r 4 /105. Since the radial functions (2.54) are orthonormal 
using 

''drrVi(r)^(r)= ( 5 nri ,, (2.61) 



one can also easily perform the integrals over r. 

2.5. Hamiltonian matrix and the eigenf unctions 

In the last section we have shown the way to construct the Hamiltonian 
matrix of the SYMQM model for D = 2 and the SU(2) group. Now, we 
are ready to construct the matrix elements and the eigenfunctions of the 
Hamiltonian in the cut Fock space (2.48). Thus, applying (2.37) and (2.59)- 
(2.61) we calculate matrix elements of (2.45) as 

H (n',i',m'Un,l,m) = ^ ^ m '\ H ^ £> m ) ? ( 2 .62) 

where \n,£, m) = \^ n ' e ' m \r, u, v)) . Next, we solve the eigenequation 

£ ff (»7,m'),W,™) B W«») = E k v^' r ' m) , (2.63) 

(n,i,m) 
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where the eigenvalues E\. describe the spectrum of the Hamiltonian (2.10) 4 . 
Finally, we find the eigenfunctions of (2.10) in the basis (2.48) as 

\<!> k (r,u,v)) = v£ Am) \¥ n > e ' m \r,u,v)) = 

(n,t,m) 

= E 4 n ''' m) E* m) (^-)^W|e,), (2.64) 

(n,l,m) m=l 

where yji' m \u,v) are spherical harmonics satisfying (2.50) while 0f t (r) are 
defined by (2.54). The non-orthogonal vectors |e M ) contain the fermionic 
variables and their scalar product is defined by (2.34). The above formula 
is the main expression which we use to compute the wave-functions of the 
model. 

3. Results 

3.1. Eigenenergies, coexistence and B — dependence 

Solving the eigenequation (2.63) for different cut-off B we obtain the 
energy spectrum of the system with j = and np = 2 as a function of B. 
Using Alpha DEC, PC computers and optimized C++ code we were able to 
reach the cut-off values B = 60. The spectrum is shown for B < 60 in Fig. 3. 
Looking at this Figure one can notice two kinds of behaviour. The levels 
from the first group are the ones which are rapidly, possibly exponentially, 
convergent to finite energy values and they become constant for large B. 
These values represent the full non-perturbative eigenspectrum of the "un- 
cut" system. The other curves fall down slowly as 

E h {B)~±. (3.1) 

These two kinds of curves do not cross but one can see, especially for 
lower B, that the levels repel each other exchanging their own -B-behaviour. 
For higher B this repulsion is more subtle. 

There is a conjecture that the asymptotics of the wave function at large 
distances, i.e. large r, determine the convergence of our calculations being 
performed with the enlarging number of allowed quanta B [12]. Thus, local- 
ized states, whose wave-functions do not go deep into the valleys, converge 
faster with the cut-off. These states are related to the discrete spectrum of 

4 To perform the calculation the van Baal's program [15] was rewritten from Mathe- 
matica code to CH — h This speeded up the program and made possible computation 
for higher cut-offs. 
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Figure 3. The energy spectrum as a function of cut-off B > ub 



60 B 



the Hamiltonian (2.10). They are supersymmetric partners of the the local- 
ized states from np = 1,3 sectors from supersymmetric multiplets [12, 19]. 
To obtain spectrum of these states the limit 



E = lim E B ,k 



fc=const 



(3.2) 



is performed where index k enumerates consecutive energy curves and changes 
only at the anti-crossing points. 

On the other hand, the non-localized states penetrate the valleys with 
increasing the cut-off. They have the power-like behaviour of the energy 
level. The energy curves which fall slowly to the zero energy with increasing 
cut-off never converge. They form the continuous spectrum at B = oo. It 
was in Ref. [18] shown that non-trivial and correct continuum limit for these 
states is given by 



E ( p ) = J im E B,k(B,P) 



with k(B,P) = -V2B (3.3) 
p=const x ' it ' 



where P is the continuous momentum, B is the cut-off and k enumerates 
only energy levels from the continuous spectrum. The continuous spectrum 
consists of all positive energy values. Localized and non-localized states 
coexist as a consequence of the supersymmetric interactions with flat valleys. 
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It turns out that the vacuum state belongs to the sector with j = and 
uf = 2 and is formed by the continuous spectrum [15, 12]. In order to obtain 
this state one has to take the states from the continuous spectrum, perform 
the continuous limit (3.3) and at the end go with E — > 0. Since the SUSY 
vacuum is constructed from the states with continuous energy spectrum, it 
is non-normalizable 5 [12]. 



4 - 




k 2 /B ~ P 2 



0.2 0.4 0.6 0.8 1. 

Figure 4. Dispersion relation, E ~ P 2 , for a few first energy curves from the contin- 
uous spectrum. Here, k - enumerates consecutive energy curves for the continuous 
spectrum and according to [18] k is related to momentum P as k/VB ~ P. 



3.2. Scaling and dispersion relation 

It was shown in Ref. [18] that for the free Hamiltonian the momentum 
of its eigenstate P is related to a discrete index of consecutive continuous- 
energy levels k by 

P ~7B' (3 " 4) 

where B is a number of bosons. It turns out that similar dispersion relation 
holds for our system. This is seen in Fig. 4 where E(P 2 ) is plotted. In this 
Figure the energies from a distinct energy curves of continuous spectrum, 



5 The vacuum state is non-normalizable means its normalization to the Dirac delta 
function. 
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enumerated by k, are denoted by different symbols. In order to evaluate P 
Eq. (3.4) was used. Indeed, the points in the plot form a curve E{P) ~ 
P 2 with a very good approximation. Thus, one can suppose that these 
states form the continuous spectrum at B = oo. This again agree with the 
continuous nature of scattering states. Therefore, we confirm that the states 
which correspond to the energy levels with power-like behaviour describe 
asymptotically free particles. One can imagine that these particles propagate 
in flat valleys of the potential (2.20). 

The above scaling law (3.3) is required to recover the infinite Hilbert 
space limit (3.3). On the other hand it is interesting because its universality. 
It is valid for models with only continuous spectrum as well as the mixed 
one, like the one considered in this work where the Hamiltonians are less 
trivial and the localized and non-localized states coexist at the same energy. 

3.3. Localized versus non-localized eigenfunctions 

Applying (2.64) one can also calculate the eigenstates of (2.45). A simple 
way to analyse the bosonic probability density is to use 

|$ fc (xi,x 2 ,x 3 )| 2 = (& k (r,u,v)\$ k (r,u,v)) F , (3.5) 

where we apply (2.34) while the integration in this product is only performed 
over the fermionic degrees of freedom. 

The energy levels of localized states converge rapidly with B. Similarly, 
their eigenstates become practically B independent. As one can see due to 
anti-crossing energy levels (see Fig. 3) the localized states with a specified 
energy (but considered for different cut-off B) corresponds to distinct values 
of k. For our convenience we denote the localized states by index (L), 

i.e. &i L \r,u,v), where contrary to & k (r, u, v) case, index s enumerates only 
states of the discrete spectrum. A few first bound states are shown in Figs. 5- 
8. There are mainly localized around the origin where potential has wide 
minimum - the stadium. These Figures show plots with various cut-off's B. 
The convergence of the first two states is so good that it is not possible to 
distinguish the wave-functions with different cut-off's. For higher states with 
s = 3, and s = 4 this convergence is worse 6 . This is caused by interference 
of the states which have nearly the same energy (see Fig. 3). The next state 
with s = 5 is not degenerate and its wave- function also rapidly converges. 

In order to observe the symmetry of the localized states we present also 
contour plot for \^ s L \xi, X2, x^)] 2 = 0.003 in Figs. 9-10. One can see that 



Because of the simpler shape, the wave-function with s = 4 converges more rapidly 
than one with s = 3. However, this convergence is much slower than one for the other 
wave-functions, i.e. with s = 1, 2, 5, . . . . 
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Figure 5. The first bound eigenstate s = 1 for B = 37, . . . , 46: 

^(z,^ 0,^ = 0) | 2 
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Figure 6. This is Fig. 5 but the bound eigenstate with s = 2 for B = 37, . . . 

\$ X s L 2 2 (x,x 2 = 0,xs = 0)\ 2 
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Figure 8. And for the fourth bound eigenstate s = 4 for B — 40, . . . , 46: 
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Figure 9. Contour plots for the bound eigenstates with s = 1,2 for 

\<fr[ L \xi,x2,xz)\ 2 = 0.0003 

only the first bound state possesses the symmetry of the bosonic potential 
(2.20): 

Xi <-» ±Xj for i, j = 1, 2, 3. (3.6) 

For the remaining states one part of the symmetry (3.6), i.e. (2.21), can 
be broken by the fermionic term in the Hamiltonian (2.33) 7 . Thus, in the 
fermionic case the symmetry (2.21) is not valid. Although in the considered 
case j = the eigenstates are not spherically symmetric and they have 
complicated shapes of contour plots with a genus sometimes greater than 
zero. However, one has to remember that the angular momentum j = 
applies for the 9— dimensional space of cf . Since the transformation (2.16) 
of the c— coordinate space to the 3— dimensional Xj-space is non-linear the 
three-dimensional spherical symmetry in the Xj-space is not obvious. 

Let us consider the probability density of states with continuous spec- 
trum along the potential valley: 

x\ = x, X2 = x% = 0. (3.7) 

It turns out that going away from x\ = X2 = X3 = the amplitude of the 
eigenstates (3.5) is strongly suppressed. Thus, to see its dependence along 



7 Hf contains v variables without square powers. This breaks the symmetry (2.21) of 
the bosonic potential. 
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Figure 10. The same contour plots but for the bound eigenstates with s = 3, 4 for 
\$r ) (x u x2,x 3 ) | 2 =0.0003 



0.0016 




-10 -5 5 ^10 



Figure 11. The functions £ 2 (x)\ 2 for eigenstates with energy levels Ek{B) 
from the continuous spectrum. Here E — 1.7 and k = 1,2, 3,4, 5. 

the valley a quantity 

x 2 |$f L \x)\ 2 = x 2 \<f> { k N L \ Xl = x,x 2 = 0,x 3 = 0)| 2 (3.8) 
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ri . TT s 0.002 

x 2 \*[ NL \x)f 

0.0018 




Figure 12. The plot similar to ones from Fig. 11 but for the eigenstates with E = 2.6 
and k = 1,2,3,4,5. 
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Figure 13. The same function as in Fig. 11 but for E = 3.7 and k = 1, 2, 3, 4, 5. 
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is studied. Similarly to the localized states the non-localized ones are de- 
noted by superscript (NL), i.e. ^^ L \xi,X2,x^), and their index k enu- 
merates only the energy levels from the continuous spectrum. The scaling 
limit implied in (3.3) is performed as follows. Take an arbitrary energy, i.e 
E = 1.7, and for a given energy curve from the continuous spectrum, i.e. 
E^ NL \b) enumerated by index k, find the cut-off B which best satisfies the 
following condition 

E[ NL) {B) = E, (3.9) 

where E is the arbitrary energy. The probability density of the Hamiltonian 
eigenstates multiplied by x 2 (3.8) is shown for E = 1.7 and k = 1,2, 3, 4, 5 in 

Fig. 11. As one can see for k = 1 the function x 2 \<&^ L \x)\ 2 has two peaks. 
For the next energy curve, i.e. k = 2, the plot has four peaks with lower 
amplitudes. Thus for the k— th curve one can see 2k peaks. In the limit 
(3.3) the amplitudes of the peaks for a given energy are nearly the same. 
This does not occur for the peaks in the centre of the potential and the 
peaks deeper sunk in the valley. For all k the oscillations of frequencies are 
nearly the same and correspond to the energy E = 1.7. For higher energies, 
i.e. E = 2.6 and E = 3.7 presented in Fig. 12 and Fig. 13 respectively, 
frequencies are lower but for the fixed energy value they are approximately 
independent of k. Taking states for higher energy curves one can notice that 
their wave- functions enter deeper into the valleys (2.20). This fact confirms 
that the considered states form the continuous spectrum for B — > oo. 

In the purely bosonic sector, i.e. where tif = 0, 6, as well as in the sectors 
with np = 1,5, where bosonic modes are dominating, despite the flatness 
of the potential valleys the flat directions are blocked by the energies of the 
transverse quantum fluctuations. This makes the spectrum in zero-fermion 
and one-fermion sectors discrete. However, in the supersymmetric system 
the transverse fluctuations cancel among bosonic and fermionic states. Thus, 
valleys are not blocked and in the model with supersymmetric fermions the 
continuous spectrum appears. This occurs in sectors for np = 2,3,4. 

On the other hand, we have sectors with only discrete spectra, i.e. 
tif = 0,1,5,6. The supersymmetry creates super-multiplets and demands 
the existence of similar states in others sectors. Therefore in the central sec- 
tors, i.e. riF = 2, 3, 4, there is coexistence of discrete and continuous spectra. 
This is exactly what Fig. 3 shows. 

4. Summary 

Despite the fact that Supersymmetric Yang-Mills Quantum Mechanics 
(SYMQM) models [1, 2] are simpler than the original field theories, they still 
pose difficult to solve complex problems. On the other hand we have the 
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BFSS equivalence hypothesis [3, 4, 5], which relates SYMQM to M— theory. 
This makes SYMQM even more interesting. Moreover, the SYMQM models 
form very good laboratory for learning of supersymmetric theories. Studies 
in their zero-fermion sectors, i.e. of the non-supersymmetric Yang-Mills me- 
chanics, give us knowledge about 0- volume glueball states, [9, 10, 11, 12]. On 
the other hand testing relations between fermionic sectors we can observe 
specific properties of the action of the supersymmetric operator. 

In this work we have considered the D = 4 dimensional model with 
SU(2) gauge groups. This model is non-trivial and possesses both localized 
and non-localized states. There were many approaches [12, 13, 15, 19] to 
understand the model and to find its energy spectrum. We have followed 
the method [15] proposed by van Baal which allows to solve problem in the 
sector with specified not only the number of fermionic quanta uf but also 
at the fixed momentum j = 0. We have chosen the sector with tif = 2 and 
j = 0, where the supersymmetric vacuum should appear. 

First, the energy spectrum for very high cut-offs B < 60 has been cal- 
culated. It exhibits complex structure of the states, i.e. coexistence of lo- 
calized and non-localized states. For a high cut-offs, i.e. B > 40, not only 
the energies of localized states become B— independent but also the corre- 
sponding eigenfunctions. We have constructed the first few bound state and 
described their properties. Our spectrum agree with the previous results 
[12, 13, 15, 19]. However, with the higher cut-off it is much more accurate. 

In the plot of the energy spectrum as a function of the cut-off B in 
Fig. 3 the non-localized states form curves behaving as 1/B. Taking the 
states from these curves for a constant energy we have studied the proper- 
ties of the eigenfunction from the continuous spectrum. At the end we have 
tested the scaling of the continuous spectrum (3.3) and confirmed numer- 
ically dispersion relation of the free particles presented in Ref. [11]. The 
scaling was derived in Ref. [18] for models with only non-localized states 
whereas here we have shown that with a very good approximation is also 
valid for the non-trivial interactions. 
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